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ABSTRACT 



o 

^""^ ' Context. In our previous work we found a statistically significant offset AV ~ 27 m s _1 between the radial velocities 

of the HC3N J = 2 — 1 and NH3 (J, K) — (1, 1) transitions observed in molecular cores from the Milky Way. This may 
O |i indicate that the electron-to-proton mass ratio, /1 = m c /m p , increases by ~ 3 x 10 -8 when measured under interstellar 

conditions with matter densities of more than 10 orders of magnitude lower as compared with laboratory (terrestrial) 
environments. 

Aims. We map four molecular cores L1498, L1512, L1517, and L1400K selected from our previous sample in order to 
estimate systematic effects in AV due to possible velocity gradients or other sources across the cloud, and to check the 
reproducibility of the velocity offsets on the year-to-year time base. 

Methods. We use the ammonia method, which involves observations of inversion lines of NH3 complemented by 
rotational lines of other molecular species and allows us to test changes in fi due to a higher sensitivity of the inversion 
frequencies to the ^-variation as compared with the rotational frequencies. 

Results. We find that in two cores L1498 and L1512 the NH3 (1,1) and HC3N (2-1) transitions closely trace the same 
material and show an offset of AV = Vi ST .(HC3N) - Vi sr (NH3) = 26.9 ± 1.2 sta t ± 3.0 sys m s _1 throughout the entire 
clouds. The offsets measured in L1517B and L1400K are 46.9 ± 3.3 s tat ± 3.0 sys m s _1 , and 8.5 ± 3.4 sta t ± 3.0 sys m s _1 , 
respectively, and are, probably, subject to Doppler shifts due to spatial segregation of HC3N versus NH3. We also 
determine frequency shifts caused by external electric and magnetic fields, and by the cosmic black body radiation- 
induced Stark effect and find that they are less than 1ms . 

Conclusions. The measured velocity offset in L1498 and L1512, being expressed in terms of Afx/fj, = (/Uobs — /Ulab)/Ati a b, 
gives Afx/fx = (26 ± l s tat ± 3 sys ) x 10 -9 . Although this estimate is based on a limited number of sources and molecular 
pairs used in the ammonia method, it demonstrates a high sensitivity of radio observations in testing fundamental 
physics. The non-zero signal in Afi/fi should be further examined as larger and more accurate data sets become 
available. 
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1. Introduction bearing molecules arising in cold and dense molecular 

clouds (T kin ~ 10 K, n ~ 10 4 - 10 5 cm" 3 ). This study is 

The present paper continues our study (Levshakov et al. aimedj in general, at testing the Einstein equivalence prin- 

2008; Molaro et al. 2009; Levshakov et al. 2010a, hereafter ciple of local posi tion invariance (LPI) which states that 

Paper I) on differential measurements of the electron-to- outcomes of physical non-gravitational experiments should 

proton mass ratio, fi = m c /m p , by means of high reso- be independent of their position in space-time (e.g., Dent 

lution spectral observations (FWHM - 30 - 40 m s x ) 2008). However, the violation of LPI is anticipated in non- 

of narrow emission lines (FWHM < 200 m s" 1 ) of N- standard physical theories associated with varying funda- 

mental constants and, in particular, in those dealing with 

* Based on observations obtained with the 100-m telescope dark energy. A concept of dark energy with negative pres- 

at Effelsberg/Germany which is operated by the Max-Planck sure appeared in physics long before the discovery of the 

Institut fur Radioastronomie on behalf of the Max-Planck- accelerating universe through observations of nearby and 
Gesellschaft (MPG). 
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distant supernovae type la (Perlmuttcr ct al. 1998; Riess 
et al. 1998). Early examples of dark energy in a form of a 
scalar field with a self-interaction potential can be found in 
a review by Peebles & Ratra (2003). To explain the nature 
of dark energy, many sophisticated models have been sug- 
gested (see, e.g., Caldwell et al. 1998), and among them the 
scalar fields which are ultra-light in cosmic vacuum but pos- 
sess an effectively large mass locally when they are coupled 
to ordinary matter by the so-called chameleon mechanism 
(Khoury & Weltman 2004a,b; Brax et al. 2004; Avelino 
2008; Burrage et al. 2009; Davis et al. 2009; Brax 2009; 
Upadhye et al. 2010; Brax & Zioutas 2010). A subclass of 
such models considered by Olive & Pospelov (2008) pre- 
dicts that fundamental physical quantities such as elemen- 
tary particle masses and low-energy coupling constants may 
also depend on the local matter density. 

In the Standard Model of particle physics there are 
two fundamental parameters with the dimension of mass - 
the Higgs vacuum expectation value (VEV ~ 200 GeV), 
which determines the electroweak unification scale, and the 
quantum chromodynamics (QSD) scale Aqcd ~ 220 MeV, 
which characterizes the strength of the strong interaction. 
The masses of the elementary particles such as the electron 
m c and the quarks m q are proportional to the Higgs VEV, 
and the mass of the composite particle proton m p is pro- 
portional to Aqcd- Therefore, by probing the dimensionless 
mass ratio \l = m e /m p we can experimentally test of the 
ratio of the electroweak scale to the strong scale at different 
physical conditions (Chin et al. 2009). 

In this regard, the masses of the elementary particles 
mediated by the chameleon-like scalar fields should depend 
on the local density of baryons and electrons. Among the 
observationally accessible media the most dramatic differ- 
ence of local densities occurs between terrestrial laborato- 
ries and cosmic space. For instance, a typical interstellar 
molecular cloud with the kinetic temperature X^m ^ 10 K 
and the gas number density n H2 ~ 10 5 cm -3 has the local 
density of p spaC c ~ 3 x 10~ 19 g cm~ 3 . Laboratory vacuum 
spectrographs usually operate in the pressure range from 1 
to 10 mtorr and at the temperature regime of T^ n ~ 10 — 40 
K (e.g., Willey et al. 2002; Ronningen & De Lucia 2005; 
Ross & Willey 2005). This gives pi ab « 3 x 10~ 9 g cm" 3 
if ?kin = 10 K, P = 1 mtorr, and the vacuum chamber is 
filled by H 2 . In this case, the ratio piab/Pspace exceeds 10 10 . 
Besides, the size of the vacuum chamber is ~ 10 cm, and 
its walls are naturally made of extremely heavy materials, 
i.e. the chameleon fields in the chamber are subject to high 
density boundary conditions. 

It follows from the above considerations that the ratio 
/i = m /m p measured in low density environments of the 
interstellar and/or intergalactic medium may differ from 
its terrestrial (laboratory) value. To estimate the fractional 
changes in fi, A/x/fi = (fi ohs - Mlab)/Mlab, the so-called 
ammonia method was elaborated by Flambaum & Kozlov 
(2007, hereafter FK) which involves observations of inver- 
sion lines of NH3 complemented by rotational lines of other 
molecular species. The inversion transitions of ammonia 
and its isotopologues strongly depend on fi (van Veldhoven 
et al. 2004; FK; Kozlov et al. 2010). For NH 3 , the fractional 
change in frequency is Afi nv /vi nv = 4.46A/z//z, whereas for 
rotational lines, Av I0 t / Vrot = A/i//i. Therefore, the com- 
parison of radial velocities of NH3 inversion lines, Vj nv , with 



radial velocities of rotational transitions, V ro ti provides a 
sensitive limit on the variation of /1 (FK): 

A/V/z « 0.289(F rot - V im )/c « 0.3AV/c , (1) 

where c is the speed of light. Thus, measuring the line shifts 
in cosmic objects where the local matter density is signifi- 
cantly lower than in laboratory experiments one can probe 
whether the mass-ratio [i is position- invariant or not. Such 
tests will complement current cosmological probes of dark 
energy and its couplings to matter. 

NH3 was detected in numerous molecular clouds within 
the Milky Way (e.g., Jijina et al. 1999) and even at sig- 
nificant redshifts (Henkel et al. 2005, 2008). In particular, 
observations of two molecular absorption-line systems at 
redshift z — 0.68 and z — 0.89 towards, respectively, the 
gravitationally lensed quasars B0218+357 and PKS 1830- 
211 yielded constraints on A/i//i down to 1 ppm (FK; 
Murphy et al. 2008; Menten et al. 2008; Henkel et al. 20090 
However, Afx/n can be probed at orders of magnitude more 
sensitive levels if nearby Galactic clouds are studied. Then 
a spectral resolution of FWHM ~ 30 m s _1 allows us 
to measure the line position with uncertainties of a few 
m s _1 . In 2008-2009, we observed 41 molecular cores in the 
disk of the Milky Way in three molecular transitions NH3 
(J,K) = (1, 1), HC 3 N J — 2-1, and N 2 H+ J = 1-0 at the 
radio telescopes in Medicina (32-m), Nobeyama (45-m), and 
Effelsberg (100-m) and obtained a statistically significant 
positive velocity offset between the rotational HC3N (2-1), 
N 2 H + J = 1 — and inversion NH3 (1,1) lines. Since these 
transitions show a good correlation in their spatial distri- 
butions, such an offset might indicate that at low interstel- 
lar densities fj, increases compared with its terrestrial value 
(Paper I). The most accurate estimate was obtained from 
the Effelsberg dataset consisting of twelve pairs of NH3 and 
HC3N lines: AV = 23 ± 4 sta t ± 3 sys m s _1 which translates 
into Afx/fj, — 22 ± 4 sta t ± 3 sys ppb. We note that because of 
a rounded value for the fast Fourier transform spectrometer 
(FFTS) channel separation (Av = 1.220 kHz for NH 3 and 
Av = 1.221 kHz for HC 3 N instead of 1.220703125 kHz) not 
taken into account in Paper I, the corrected A[i/[i is to be 
26 ± 4 stat ± 3 sys ppb. 

These observations were carried out in a single position 
mode, i.e. we observed only the central parts of molecu- 
lar cores showing strongest NH3 emission. The aim was to 
find suitable targets for most precise measurements of the 
velocity offsets AV^ between rotational and inversion tran- 
sitions. In the present paper, we report on new observa- 
tions where we measure AV at different positions across 
individual clouds in order to test the reproducibility of 
the velocity offsets in the presence of the large-scale veloc- 
ity gradients. For this purpose, from the list of molecular 
cores observed in our previous study with the Effelsberg 
telescope, we selected the objects which comply with the 
following criteria: (i) profiles of the NH3 (1,1) and HC3N 
(2-1) hyperfine structure (hfs) transitions are symmetric 
and well-described by a single-component Gaussian model, 
and (ii) the line widths are thermally dominated, i.e., the 
parameter /3 = f7 t ,(NH 3 )/(T„(HC3N) > 1, where a v is the 
velocity dispersion. The chosen targets are molecular cores 
L1498, L1512, L1517B and L1400K which were already ex- 
tensively studied in many molecular lines (e.g. Benson & 
Myers 1989; Lee et al. 2001; Lee et al. 2003; Tafalla et al. 

1 Hereafter, 1 ppm = 1CT 6 , and 1 ppb = 10~ 9 . 
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Table 1. Source positions, offsets (Aa,A6), and date of observations. 
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2004, 2006; Crapsi et al. 2005). In the present paper we 
describe the high-precision measurements of the radial ve- 
locities of NH3 (1,1) and HC3N (2-1) hyperfme transitions 
in these clouds from 31 lines of sight in total. 



2. Observations 

We carried out our observations targeting a sample of well- 
studied nearby molecular cores that are essentially devoid of 
associated IR sources (so-called starless cores with narrow 
emission lines) . All of our objects lie in the Taurus- Auriga- 
Perseus molecular complex, i.e., at the distance of about 
140 pc (L1498, L1512, L1517BC) and 170 pc (L1400K) 
from the solar system (Jijina et al. 1999). Observations 
were performed with the 100-m Effelsberg radio telescope in 
January 2010. The positions observed are listed in Tabled] 
All sources, except L1517B, have the same central coor- 
dinates as used in Paper I. For L1517B, the central posi- 
tion is shifted at (Aa, AS) — (15", 15") with respect to the 
previous coordinates, and the source is called L1517BC. 
The (J,K) — (1,1) inversion line of ammonia (NH3) at 
23.694 GHz and the J = 2 — 1 rotation line of cyanoacety- 
lene (HC3N) at 18.196 GHz were measured with a K- 
band HEMT (high electron mobility transistor) dual chan- 
nel received, yielding spectra with angular resolution of 
HPBW ~ 40" in two orthogonally oriented linear polar- 
izations. Averaging the emission from both channels gives 
typical system temperatures of 100-150 K for NH3 and 80- 
100 K for HC3N on a main beam brightness temperature 
scale. 



2 At 18 GHz we used only one channel Rx2 since Rxl was not 
operating during our observations. 



The measurements were carried out in frequency switch- 
ing mode using a frequency throw of 5 MHz. The backend 
was an FFTS operated with its minimum bandwidth of 
20 MHz providing simultaneously 16 384 channels for each 
polarization. The resulting channel separations are 15.4 
m s _1 and 20.1 m s _1 for NH3 and HC3N, respectively. 
We note, however, that the true velocity resolution is about 
twice as lower, FWHM ~ 30 m s" 1 and 40 m s , respec- 
tively (Kein et al. 2006). 

The sky frequencies were reset at the onset of each scan 
and the Doppler tracking was used continuously to track 
Dopplcr shifts during the observations. In the Effelsberg 
driving program, the position of the center of gravity of 
Earth is taken from the JPL ephemcrides including the in- 
fluence of the Moon, and the rotation of Earth at the tele- 
scope position is added to this center of gravity. The ve- 
locity of the source with respect to Local Standard of Rest 
(LSR), Visr, is corrected for these motions and the corrected 
value of Vi sr determines the current sky frequency. The sta- 
bility of the LO-frcqucncy during 5 minutes exposures was 
about 10" 14 . 

Observations started by measuring the continuum emis- 
sion of calibration sources (NGC 7027, W3(OH), 3C286) 
and continued by performing pointing measurements to- 
ward a source close to the spectroscopic target. Spectral 
line measurements were interspersed with pointing mea- 
surements at least once per hour. The calibration is es- 
timated to be accurate to ±15% and the pointing accu- 
racy to be superior to 10 arcsec. Each molecular transition 
was observed in a set of four sequential 5 min exposures, 
which were repeated twice at 18 GHz to compensate for 
the non-operating channel Rxl. Thus, the total exposure 
time at one coordinate point was one hour. Some points 
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Fig. 1. Maps of L1498, L1512, L1517BC, and L1400K showing positions where the radial velocity differences, AV, 
between the HC 3 N (2-1) 18.196 GHz line and the NH 3 (1,1) 23.694 GHz line have been measured at the Effelsberg 
100-m radio telescope. On each panel the numbers are the mean velocity shifts AV (in m s _1 ) between the corresponding 
AV S and AV a values listed in Tables [5HS] The circles show the half-power beam width at 23 GHz, whereas the linear 
scale is indicated by the horizontal bars. Reference positions are given in Table [1] 



were observed several times to control the accuracy of the 
AV measurements. All spectra were obtained with a com- 
parable noise level which is indicated by the rms values in 
Table [T] The position offsets within each individual cloud 
are shown in Fig. [T] 

3. Analysis 

We used the CLASS reduction packageQ for standard data 
reduction. After corrections for the rounded frequencies (see 
Sec.rjJ, the individual exposures were co-added to increase 
the signal-to- noise ratio, S/N. The spectra were folded to 



remove the effects of the frequency switch, and base lines 
were determined for each spectrum. The resolved hfs com- 
ponents show no kinematic sub-structure and consist of an 
apparently symmetric peak profile without broadened line 
wings or self-absorption features, as shown in Fig. [2] where 
the observed profiles of the NH 3 (1,1) and HC 3 N (2-1) lines 
toward L1498 are plotted. The line parameters such as the 
total optical depth in the transition, r to t (i.e., the peak op- 
tical depth if all hyperfine components were placed at the 
same velocity), the radial velocity, Vi sr , the line broadening 
Doppler parameter, b, and the amplitude, A, were obtained 
through the fitting of the one-component Gaussian model 
to the observed spectra as described in Paper I: 



3 http://www.iram.fr/IRAMFR/GILDAS 



T(v) = A- [1 - exp(-t(v))} , 



(2) 
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Table 2. Molecular core L1498. Represented arc fitting parameters Vi sr (km s -1 )/b (km s -1 )/ xt- Lower limits on Tkin, 
and upper limits on <7t U rb are indicated. The numbers in parentheses are la statistical errors. The signal-to-nose ratio, 
S/N, per spectral channel at the maximum intensity peak of NH3/HC3N is indicated in the last column. 
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Note. — Numbering in column 2 is the same as in the last column of Table [T] 



with 



= T ■ X! ai eXP i~( V ~ V i~ V lsr) 2 /b 2 ] , 



which transforms for optically thin transitions into 



T(v) = A' 



exp [-(v -Vi — Vi sr ) 2 /b 2 



(3) 



(4) 



The sum in and (H|) runs over the k = 18 and k = 6 hfs 
components of the NH 3 (J, A") = (1, 1) and HC 3 N J = 2 - 1 
transitions, respectively. 

The rest-frame frequencies of the NH3 hfs components 
are taken from Kukolich (1967) and listed in Table 2 
(Paper I). We checked these frequencies using the last JPL 
modcQ which takes into account the fit to all published am- 
monia data. The difference between various estimates does 
not exceed 0.05 kHz. In addition, Lapinov et al. (2010) car- 
ried out new laboratory measurements of the frequencies 
of the NH3 (1,1) hfs components and their results coincide 
with Kukolich's frequencies within the 0.2 kHz range. 

The rest-frame frequencies of the HC3N (2-1) hfs com- 
ponents listed in Table 3 (Paper I) are based on Lapinov 
(2008, private comm.) calculations which combine all avail- 
able laboratory measurements. These frequencies are in 
agreement with the data from the Cologne Database for 
Molecular Spectroscopy (Muller et al. 2005) within la- 
uncertainty interval, but Lapinov's values have a slightly 
higher precision: e„ ~ 0.17 kHz versus 0.2 to 0.7 kHz by 
Muller et al.. The recent laboratory measurements of the 
HC3N (2-1) frequencies by Lapinov et al. (2010) did not 



show any systematic shifts between the previously obtained 
and new measured values: all frequencies agree within 0.2 
kHz. Note that at the moment the accuracy with which 
the frequencies of NH3 (1,1) and HC3N (2-1) are known 
is the highest among all observable molecular transitions, 
with the uncertainties being £„(NH3) ~ 0.6 m s _1 (~ 50 
Hz), and e„(HC 3 N) ~ 2.8 m s^ 1 (~ 170 Hz). 

We also estimated the kinetic temperature, Tk; n , and 
the nonthermal (turbulent) velocity dispersion, <7t U rb, based 
on the line broadening Doppler parameters, b — \/2a v , of 
the NH 3 (1,1) and HC 3 N (2-1) lines. Here a v is the line 
of sight velocity dispersion of the molecular gas within a 
given cloud. If the two molecular transitions trace the same 
material and have the same non-thermal velocity compo- 
nent, then a v is the quadrature sum of the thermal <7th arid 
turbulent dturb velocity dispersions. In this case a lighter 
molecule with a mass mi should have a wider line width as 
compared with a heavier molecule with ro/, > mi. 

For thermally dominated line widths (f3 > 1, see Sect. 1) 

and co-spatially distributed species one has the following 

relations (e.g., Fuller & Myers 1993): 

mim h 2 2 

Oz -o- h ), (5) 



T, 



kin 



and 



k(m h - mi) 



2 _ m h <rl ~ miaf 

T turb 



http:/ /spec.jpl. nasa.gov/ftp/pub/catalog/catdir. html 



(6) 

m h - mi 

where k is Boltzmann's constant, and the thermal velocity 
dispersion, cr t h, is given by 

<7th,i = (fcTkin/m,) 172 . (7) 

However, HC3N is usually distributed in a larger volume of 
the molecular core as compared with NH3: in general, N- 
bearing molecules trace the inner core, whereas C-bearing 



6 S. A. Levshakov et al.: Searching for chameleon-like scalar fields with the ammonia method 

Table 3. Same as Table [5] but for the molecular core L1512. 







JNH3 


(I. 1 ) 


IT/" TVT /O -I \ 

HC3JN (2—1) 


AT/ 

A 14, 


AT/ 


J kin, 


Oturb, 


C /T\T 
B/lN 


Offset 


# 


satellite 


all 


all 


m s _1 


m s _1 


K 


m s _1 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0,0 


1 


7.1015(18) 

0.112(5) 

1.17 


7.1015(14) 

0.110(3) 

1.11 


7.1290(21) 

0.075(9) 

0.94 


27.5(2.8) 


27.5(2.5) 


10.0(2.3) 


34.6(14.8) 


24 
19 


0,0 


10 


7.1035(12) 

0.108(4) 

1.13 


7.1045(9) 

0.111(1) 

1.11 


7.1300(20) 

0.075(5) 

1.13 


26.5(2.3) 


25.5(2.2) 


10.3(1.2) 


33.7(8.4) 


35 
20 


0,0 


12 


7.1060(12) 

0.106(3) 

1.08 


7.1070(9) 

0.111(1) 

1.17 


7.1320(18) 

0.081(7) 

1.16 


26.0(2.2) 


25.0(2.0) 


8.9(1.8) 


42.9(9.9) 


27 
16 


0,0 


13 


7.1050(13) 

0.109(4) 

0.95 

7.1030(85) 

0.113(12) 

0.94 

7.0380(9) 

0.110(3) 

1.07 


7.1070(12) 

0.111(3) 

0.93 

7.1150(52) 

0.112(12) 

0.95 

7.0385(8) 

0.109(1) 

1.16 


7.1295(20) 

0.071(10) 

1.13 

7.1475(52) 

0.071(8) 

1.07 

7.0710(29) 

0.091(22) 

1.06 


24.5(2.4) 


22.5(2.3) 


11.2(2.4) 


26.5(20.4) 


29 
16 


20,-40 


2 


44.5(10.0) 


32.5(7.4) 


11.6(4.5) 


25.4(21.4) 


9 
7 


-20,40 


3 


33.0(3.0) 


32.5(3.0) 


5.6(6.2) 


56.9(26.4) 


38 
15 


-40,80 


4 


7.0130(22) 

0.116(5) 

0.92 


7.0115(17) 

0.117(3) 

0.93 


7.0230(43) 

0.084(13) 

0.70 


10.0(4.8) 


11.5(4.6) 


10.2(3.5) 


43.2(19.0) 


22 
10 


-30,60 


5 


7.0225(14) 

0.116(4) 

1.09 


7.0215(10) 

0.116(1) 

1.17 


7.0475(48) 

0.090(38) 

1.02 


25.0(5.0) 


26.0(4.9) 


8.3(10.6) 


52.1(49.3) 


32 
14 


-10,20 


6 


7.0740(10) 

0.110(5) 

0.91 


7.0755(8) 

0.107(1) 

1.07 


7.1075(23) 

0.090(5) 

0.87 


33.5(2.5) 


32.0(2.4) 


5.2(1.4) 


56.7(6.0) 


38 
18 


10,-20 


7 


7.1200(25) 

0.104(7) 

1.09 


7.1225(18) 

0.106(4) 

1.11 


7.1370(23) 

0.076(5) 

1.09 


17.0(3.4) 


14.5(2.9) 


8.4(1.8) 


39.0(7.8) 


19 
15 


-40,0 


8 


7.0935(12) 

0.107(3) 

1.01 


7.0960(9) 

0.110(1) 

1.16 


7.1300(21) 

0.098(4) 

1.05 


36.5(2.4) 


34.0(2.3) 


3.8(1.3) 


64.6(4.6) 


33 
21 


0,40 


9 


7.0505(14) 

0.105(4) 

0.76 


7.0500(12) 

0.108(3) 

0.90 


7.0860(23) 

0.080(5) 

1.10 


35.5(2.7) 


36.0(2.6) 


8.1(1.6) 


43.4(7.2) 


27 
19 


-50,100 


11 


6.9875(78) 

0.131(20) 

0.91 


6.9720(49) 

0.141(7) 

1.00 


6.9985(53) 

0.088(10) 

0.70 


11.0(9.4) 


26.5(7.2) 


18.7(4.1) 


28.9(24.3) 


11 

7 



molecules occupy the outer part (e.g., Di Francesco et al. 
2007). Such a chemical differentiation and velocity gradi- 
ents within the core may cause a larger nonthermal compo- 
nent in the velocity distribution of HC3N. If both molecules 
are well shielded from the external incident radiation and 
the gas temperature is mainly due to the heating by cosmic 
rays, then a formal application of Eqs.JS]) and © to the 
apparent line widths provides a lower limit on T^ n and an 
upper limit on cr tur b. 

We now consider the relative radial velocities of HC3N 
and NH 3 , = V rot (HC 3 N) - T4 nv (NH 3 ). As discussed in 
Paper I, the velocity offset AV^ can be represented by the 
sum of two components 

AV = AV IM + AV n , (8) 

where AV^ is the shift due to a putative ^-variation, and 
AV n is the Doppler noise — a random component caused 
by possible spatial segregation of HC3N versus NH3 and 
their different irregular Doppler shifts. The Doppler noise 
can mimic or obliterate a real signal and, hence, should be 
minimized. This can be achieved either by detailed map- 
ping of the velocity field in molecular lines which closely 
trace each other as performed in the present study, or by 
averaging over a large data sample as realized in Paper I. 



4. Results 

Equations (J2J)- (J4j) define our model of the hfs line profiles. 
The model parameters were determined by using a non- 
linear least-squares method (LSM). The best fit values of 
the radial velocity Vi sr and the Doppler width b are given 
in Tables dJini columns 3, 4, and 5. An example of the fit 
is shown in Fig. 2. The residuals depicted below the pro- 
files are normally distributed with zero mean and disper- 
sion equal to the noise dispersion in the observed spectrum. 
While dealing with NH3, we calculated two sets of the fit- 
ting parameters: (?) based on the analysis of only optically 
thin satellite lines with AFi 7^ (col. 3), and (ii) obtained 
from the fit to the entirety of the NH3 (1,1) spectrum in- 
cluding the main transitions with Ai 7 ! = which have opti- 
cal depths t > 1 (col. 4), as can be inferred from the relative 
intensities of the hfs components. The corresponding radial 
velocity offsets are marked by A V s and A V a ■ All hfs transi- 
tions of the HC3N (2-1) line from our sample are optically 
thin (r < 1) and their Vi sr and b values are listed in col. 5. 
The calculated limiting values of the kinetic temperature 
Ikm (lower limit) and the nonthermal rms velocity disper- 
sion CTturb (upper limit) of the material presumably traced 
by the NH3 and HC3N emission are given in columns 8 and 
9. Column 10 presents the signal-to-noise ratio per spectral 
channel at the maximum intensity peak of NH3/HC3N. 
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NH 3 


(1,1) 


HC 3 N (2 - 1) 




AVa, 


Tkin, 


""turb, 


S/N 


Offset 


# 


satellite 


all 


all 


m s _1 


m s _1 


K 


m s _1 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0,0 


1 


5.7975(9) 

0.115(3) 

1.11 


5.7990(8) 

0.118(1) 

1.09 


5.8515(26) 

0.111(9) 

0.95 


54.0(2.8) 


52.5(2.7) 


2.5(3.1) 


75.9(9.9) 


42 
18 


0,0 


6 


5.7935(9) 

0.114(3) 

1.13 


5.7945(8) 

0.116(1) 

1.09 


5.8490(21) 

0.098(7) 

1.06 


55.5(2.3) 


54.5(2.2) 


5.9(2.1) 


62.0(8.3) 


37 
23 


15,15 


2 


5.7835(9) 

0.118(3) 

1.21 


5.7830(8) 

0.118(1) 

1.10 


5.8250(31) 

0.106(10) 

0.71 


41.5(3.2) 


42.0(3.2) 


4.2(3.3) 


70.3(11.3) 


43 
15 


-15,-15 


3 


5.8075(9) 

0.115(3) 

1.04 


5.8085(8) 

0.115(1) 

1.14 


5.8500(27) 

0.105(8) 

0.86 


42.5(2.8) 


41.5(2.8) 


3.4(2.6) 


70.4(9.0) 


42 
18 



15,-15 
-15,15 



5.7930(9) 

0.116(4) 

0.92 

5.7990(13) 

0.119(5) 

0.87 



5.7930(8) 

0.114(1) 

1.11 

5.8020(10) 

0.114(1) 

1.07 



5.8325(33) 

0.107(9) 

0.88 

5.8535(25) 

0.099(9) 

1.01 



39.5(3.4) 39.5(3.4) 2.3(3.0) 73.1(9.9) 



41 
13 



54.5(2.8) 51.5(2.7) 4.9(2.8) 64.0(10.4) 32 

18 



The lcr errors of the fitting parameters were calculated 
from the diagonal elements of the covariance matrix at the 
minimum of \ 2 - In addition, the error of the radial veloc- 
ity Vi sr was calculated independently by the A% 2 method. 
However, because of the spectral leakage (Klein et al. 2006), 
the fluxes in neighboring spectral channels are correlated 
(correlation coefficient r ~ 0.61), and these errors must be 
further corrected. The correlation length of £ = 2 spec- 
tral channels can be deduced from the analysis of the au- 
tocorrelation function (ACF) of the intensity fluctuations 
in the spectral regions free from emission lines (computa- 
tional details are given in Levshakov et al. 1997, 2002). The 
data become less correlated (£ = 1, r ~ 0.36) after aver- 
aging over two neighboring spectral channels. The LSM is 
known to be an unbiased and minimum variance estimator 
(Gauss-Markov theorem) irrespective of the distribution of 
the measurements including correlated data as well (Aitkcn 
1934). This means that the best fit model parameters de- 
duced from both the strongly and less correlated spectra 
should coincide. Numerical tests confirm this assumption: 
from both spectra we obtained exactly the same Vi sr and 
b values, but in case of strongly correlated data the errors 
representing the A\ 2 — 1 ellipsoid were a factor of ~ 1.3 
smaller as compared with the errors based on less corre- 
lated data. Thus, errors given in parentheses in Tables [5][5l 
columns 3-5, represent calculated values multiplied by this 
factor. These errors were further transformed into uncer- 
tainties of the quantities listed in columns 6-9. We also 
checked both folded and unfolded spectra and found that 
the measured velocity offsets arc not affected by the folding 
procedure. 

Since HC3N was observed with one K-band channel 
(Rx2) and NH3 with both (Rxl and Rx2), we compared 
the Rx2 data of HC3N with NH3 from the same channel 
and found that changes in the radial velocities were in- 
significant, within la uncertainty range. 

Besides the formal statistical (model fitting) uncertain- 
ties there can be additional errors caused by instrumental 
imperfections. To control this type of errors, we carried out 
repeated observations at 10 offsets having the same coordi- 
nates: #3, 8 (L1498), 1, 10, 12, 13 (L1512), 1, 6 (L1517BC), 
and 1, 2 (L1400K) in Table ED The A^ dispersions re- 



sulting from these repeated measurements (Tables OS]) ar e 
a(AV s ) = 1.6 ms" 1 andcr(AV a ) = 2.0 m s" 1 . These disper- 
sions are systematically lower than the reported la errors 
of AV S and A V a , which ensures that we are not missing any 
significant instrumental errors at the level of a few m s _1 . 

The comparison of the velocity dispersions determined 
from the NH 3 (1,1) and HC 3 N (2-1) lines (Tables^© does 
not show any significant variations with position within 
each molecular core. All data are consistent with thermally 
dominated line broadening, i.e., as expected, 1 < f3 < y/3 
for both sets 's' and 'a' of the data points. In particu- 
lar, the latter set gives the following weighted mean values: 
/3 L1498 = 1.24 (3), /3 L1512 = 1.32(5), /3 L1517BC = 1.11(1), and 
/?li4ook = 1.23(4). The maximum contribution of the non- 
thermal motions to the velocity dispersion is observed in the 
L1517BC core. The weighted mean values of the velocity 
dispersions (b/y/2) for NH 3 range between a^ l5l2 — 78(1) 
m s _1 and a L1400K — 86(1) m s _1 , and for HC3N between 
c L i5i2 = 59(2) m s _1 and a L151VBC — 74(1) m s _1 . This can 
be compared with the speed of sound inside a thermally 
dominated region of a cold molecular core which is defined 
as (e.g., Shu 1977) 

v s = (fcTkin/mo) 1 / 2 , (9) 

where too is the mean molecular mass. With too ~ 2.3 
amu for molecular clouds, one has v s ~ 60VTki n m s _1 , 
which shows that at the typical kinetic temperature of 10 
K, the nonthermal velocities are in general subsonic, and 
that the selected targets do represent the quiescent material 
at different distances from the core centers. 

The nonthermal velocity dispersions derived from the 
apparent line widths depend on the adopted gas temper- 
ature. The gas temperature in a molecular cloud is deter- 
mined by the balance between heating and cooling. If the 
only source of heating is the cosmic rays and the cool- 
ing is due to the line radiation mainly from CO, then 
a lower bound on the kinetic temperature is about 8 K 
(Goldsmith & Langer 1978). In cloud cores where the gas 
is well-shielded from the background ionizing radiation the 
gas temperature is about 10 K (Goldsmith 2001). In partic- 
ular, just this value is measured in the L1498, L1517B and 
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Table 5. Same as Table [5] but for the molecular core L1400K. 







NH 3 


(1,1) 


HC 3 N (2 - 1) 


AV a , 


AV a , 


Tkin, 


""turb, 


S/N 


Offset 


# 


satellite 


all 


all 


m s _1 


m s _1 


K 


m s _1 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0,0 


1 


3.2535(27) 

0.121(g) 

1.15 


3.2540(20) 

0.122(3) 

1.17 


3.2610(40) 

0.104(16) 

0.98 


7.5(4.8) 


7.0(4.5) 


6.3(5.3) 


66.3(18.9) 


20 
12 


0,0 


2 


3.2485(25) 

0.112(7) 

0.86 


3.2475(20) 

0.119(3) 

0.85 


3.2465(34) 

0.103(7) 

0.98 


-2.0(4.2) 


-1.0(3.9) 


5.5(2.5) 


66.5(8.2) 


19 
13 


0,-40 


3 


3.2135(26) 

0.132(5) 

1.16 


3.2160(21) 

0.125(3) 

1.11 


3.2270(26) 

0.099(5) 

0.98 


13.5(3.7) 


11.0(3.3) 


9.0(1.9) 


58.7(6.5) 


18 
17 


40,40 


4 


3.3540(27) 

0.126(7) 

1.01 


3.3570(20) 

0.120(3) 

1.10 


3.3720(31) 

0.079(16) 

1.08 


18.0(4.1) 


15.0(3.7) 


12.6(4.1) 


32.9(29.0) 


22 
15 



L1512 cores considered here (Benson & Myers 1989; Lee et 
al. 2003; Tafalla et al. 2004). 

In the L1498 cloud, we obtain a lower limit on the 
kinetic temperature Tk; n = 7.1 ± 0.5 K (average over 8 
points) which is slightly lower than T^ n = 10 K mea- 
sured by a different method from the relative population 
of the (J, K) = (2, 2) and (1,1) levels of NH 3 described by 
the rotational temperature T^ 1 (Tafalla et al. 2004). For 
L1512, the temperature averaged over 11 points (Table [3]) 
is Ttin = 9-6 ± 0.6 K which is consistent with the value 
of 10 K. The kinetic temperatures in the L1517BC core is 
well below the Tafalla et al.'s value of 9.5 K in all scanned 
points which means that the non-thermal velocity disper- 
sions of NH3 and HC3N differ significantly and that both 
species do not trace the same material. The three scanned 
points in the L1400K core show Tkin = 8.3±3.2 K - close to 
the expected value (measurements of the gas temperature 
in this core were not performed in previous studies). We 
note, however, that when the linewidths of NH3 and HC3N 
are comparable, the estimate by means of Eq. (JSJ becomes 
unstable and leads to unphysical results (e.g., points #6, 
8 in Table [3]) . In general, the spatial fluctuations of Tki n 
do not exceed a few kelvin implying the uniform heating 
and the absence of the localized heat sources. The kinetic 
temperature tends to rise with the distance from the core 
center (points #7 in Tabled #H in Tabled and #4 in 
Table 0|) which is in line with the results of Tafalla et al. 
(2004). 

The radial velocity profiles along the different diago- 
nal cuts toward the selected targets are depicted in Figs. [3] 
and [U The diagonal cut in L1498 between r — —80" and 
r = 40" (Fig. exhibits coherently changing velocities 
of V/ sr (HC 3 N) and Vi sr (NH 3 ) (except the point r = 40" 
at the core edge) with a turn at r = 20". The velocity 
gradient is small, |VV/ sr | « 0.5 km s _1 pc -1 . This picture 
coincides with the previously obtained results based on ob- 
servations of CO, CS, N2H+, and NH3 in this core and was 
interpreted as an inward flow (Lee et al. 2001; Tafalla et 
al. 2004). Taking together, all available observation classify 
the L1498 cloud as one of the most quiet molecular core. 
Thus, we can expect that in this core the Doppler noise, i.e., 
the irregular random shifts in the radial velocities between 
different transitions, is minimal. 

In the L1512 core, the VJ sr (NH 3 ) and V/ sr (HC 3 N) dis- 
tributions are almost parallel (Fig.[3)D, d). The same kine- 
matic picture was obtained for this core by Lee et al. (2001) 
from observations of CS and NH3 lines and interpreted as 



a simple rotation around the center. The velocity gradients 
derived from both the NH3 and HC3N lines are similar, 
VV/ sr ~ 1.5 km s _1 pc -1 , and consistent with the gradi- 
ent based on N2H+ measurements by Caselli et al. (2002). 
This means that NH3, HC3N, and N2H + trace the same 
gas and the Doppler shifts AV^ between them should be 
insignificant. 

In the core L1517BC which is known to be very com- 
pact (Lee et al. 2001) we observed only a small central 
part 30" x 30" where the velocities of NH 3 and HC 3 N 
along two perpendicular cuts (Fig. 2k, b) do not change 
much: VV/ sr (NH 3 ) « 0.3 km s" 1 pc" 1 , VVj sr (HC 3 N) « 0.8 
km s- 1 pc- 1 (panel a), and VVL r (NH 3 ) « W /sr (HC 3 N) 
« —0.8 km s -1 pc -1 (panel b). However, a wider area 
(« 80" x 80") observation of this core revealed an outward 
gas motion at the core periphery with a higher velocity 
gradient, W /sr (N 2 H+) « 1.1 km s _1 pc~* (Tafalla et al. 
2004) which is consistent with earlier results on NH 3 ob- 
servations by Goodman et al. (1993). This can lead to an 
additional shift in the radial velocity of HC3N line since, 
in general, the C-bearing molecules occupy a lower density 
gas (n ~ 10 3 cm -3 ) in the envelope of the molecular core. A 
higher non-thermal velocity dispersion of the HC3N line as 
compared to NH3 was already mentioned above in regard 
with the temperature measurements in this core. 

In the L1400K core we observed only three positions. 
Both molecules trace the same gradient of VVi sr ~ 1.9 
km s" 1 pc" 1 which is in line with VV/ ST -(N 2 H+) = 1.8±0.1 
km s _1 pc -1 derived by Caselli et al. (2002). The map- 
ping in different molecular lines by Tafalla et al. (2002) 
revealed that L1400K deviates significantly from spherical 
symmetry and exhibits quite complex kinematic structure. 
In particular, distributions of N2H 4 " and NH3 do not coin- 
cide: N2H + has an additional component to the west from 
the center. This explains why Craspi et al. (2005) report 
for this core a blue-ward skewness, 8 — —0.42 ±0.10, of the 
the N2H + (1-0) hfs profiles, whereas in our observations 
the NH3 (1,1) hfs transitions are fully symmetric: at the 
central position, the skewness of the resolved and single hfs 
component F\F = ± 1 1± of NH 3 is 9 = -0.1 ± 0.3. 

The mean values of the velocity offsets AV between 
HC3N and NH3 measured across each of 4 cores are pre- 
sented in Table[S] For each individual molecular core we cal- 
culated the unweighted and weighted mean AV^ (weights in- 
verse proportional to the variances) by averaging over avail- 
able offsets for two datasets: based on only optically thin 
hfs transitions of NH3 (marked by index s), and on the all 
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J/iL-A-V WApJ^ ^°^4j^4^V W^hA^ 

. . I H ... I .. I .... I ... I . . . I .... I .... I .... I .... I . . . I . . . 



■13 -12 -11 -1 



15 16 



27 28 



HC 3 N (2-1) 



F = 3-2 




F=2-l 



F=l -0 



F=2-2 



-29 



5 6 7 8 9 25 

Radial Velocity, V (km s- 1 ) 




Fig. 2. Spectra of NH3 (1,1) and HC3N (2 — 1) lines toward the molecular core L1498 obtained at the Effelsberg 100-m 
radio telescope. The offset (Aa, A6) = (0",0"). The histograms show the data averaged over two spectral channels. The 
fits are shown with solid curves. The vertical tick-marks indicate the positions and the relative strengths of the hyperfine 
structure components in case of local thermodynamical equilibrium and optically thin emission. The residual spectra are 
plotted below each profile. The fitting parameters are presented in Tabled] 



transitions including the strong main components of NH3 
(marked by index a). We also used a robust redescending 
M-estimate (Maximum-likelihood) for the mean and the 
normalized Median Absolute Deviation (1.483MAD) for 
the scale of the distribution as described in the Appendix of 
Paper I. These statistics work well for inhomogeneous data 
sets with outliers and deviations from normality. The me- 
dian values, which are robust estimates independent on the 
assumed distribution of data, are also presented to demon- 



strate the consistency of different calculations of the mean 
AV^ values. 

We obtain very similar velocity shifts (AV a ) = 25.8±1.7 
m s _1 and 28.0 ± 1.8 m s _1 (M-estimates) for, respectively, 
the cores L1498 and L1512 where the minimal level of the 
Doppler noise is expected. A larger shift (AV Q ) = 46.9 ±3.3 
m s _1 is observed in the L1517BC core - again in accord 
with the revealed kinematic structure of this core which 
allows us to expect a higher radial velocity for the HC3N 
line. On the other hand, a lower value (AV a ) — 8.5 ± 3.4 
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Fig. 3. The line of sight velocities of NH3 (J, K) = (1,1) (squares) and HC3N J = 2 — 1 (triangles) at different radial 
distances along the main diagonal cuts (panels a and b) and in the perpendicular directions (panels c and d) of the 
molecular cores L1498 and L1512 shown in Fig. \T\ The reference point in panels a and b is (Aa, A6) = (0", 0"), but it is 
(20", 0") and (—10", 20") in panels c and d, respectively. The radial distances of the points from the circular sector with 
the central angle 0° < ip < 180° are positive, whereas those from the sector 180° < ip < 360° are negative. For multiple 
observations at the same coordinate the mean values of Vi sr are depicted. The error bars show the la uncertainties. 




Fig. 4. The line of sight velocities of NH3 (J,K) = (1,1) (squares) and HC3N J = 2 — 1 (triangles) at different radial 
distances along the perpendicular diagonal cuts of the molecular core L1517BC (panels a and b), and along three points 
of L1400K (panel c) shown in Fig. [TJ The origin of coordinates is (Aa, Ad) = (0", 0") in each panel. The radial distances 
of the points from the circular sector with the central angle 0° < ip < 180° are positive, whereas those from the sector 
180° < ip < 360° are negative. For multiple observations at the same coordinate the mean values of Vi sr are depicted. 
The error bars show the la uncertainties. 
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Fig. 5. Comparison of the relative radial velocities between 
the HC3N (2-1) and NH3 (1,1) lines measured at the same 
coordinates in Feb 2009 (filled circles) and in Jan 2010 
(open circles). The 2009 data are corrected by 4 m s" 1 as 
described in Sect.[Q For L1517B the offset (A«, A<5) 2 oio = 
(15", 15") corresponds to the offset (Aa, A<5) 20 09 = (0", 0"). 
For multiple observations at the same coordinate the mean 
value is depicted. The error bars show the la uncertainties. 



m s _1 in L1400K may be due to irregular kinematic struc- 
ture of the core center which could increase the radial ve- 
locity of the NH3 line. Thus, as the reference velocity offset 
we choose the most robust Af-estimate of the mean value 
from the L1498 and L1512 cores: (AV a ) = 26.9 ± 1.2 stat 
m s . Taking into account that the uncertainty of the 
HC3N (2-1) rest frequency is about 3 m s , whereas 
that of NH3 (1,1) is less than 1 m s — 1 , we finally have 
(AV a ) = 26.9 ± 1.2 sta t ± 3.0 sys m s _1 . Being interpreted 
in terms of the electron-to-proton mass ratio variation, as 
defined in Eq. (JTJ) , this velocity offset provides Afi/fi = 
26 ± l s tat ± 3 sys ppb. 

In Fig. [5l we compare the velocity offsets AV obtained 
in our observations with the 100-m Effelsberg telescope in 
Feb 2009 and Jan 2010. For L1498, L1512 and L1517BC, 
the reproducibility is very good, whereas for L1400K the 
results differ significantly. The reason for this discrepancy 
is unclear. 



5. Discussion 

5.1. Kinematic velocity shifts 

In two molecular cores, L1498 and L1512, with lowest 
Doppler noise, we register very close values of the veloc- 
ity offset AV ~ 27 m s _1 between the rotational transition 
HC3N (2-1) and inverse transition NH3 (1,1), and these val- 
ues coincide with the most accurate estimate obtained from 
the Effelsberg dataset on 12 molecular clouds in the Milky 
Way (Paper I). Two other cores, L1517B and L1400K, ex- 
hibit velocity shifts which are either higher (~ 47 m s _1 in 
Lf517B) or lower (-9 m s _1 in L1400K) than the mean 
value, but the positive (L1517B) and negative (L1400K) de- 
flections from the mean can be explained from the observed 
kinematics in these cores. 

The rotational transition HC3N (2-1) was chosen for 
several reasons: its laboratory frequency is known with an 
accuracy of — 3 m s -1 , it closely traces the NH 3 (1,1) 
emission, and its frequency of 18.2 GHz is close to 23.7 
GHz of NH3 (1,1) yielding the similar angular resolution for 
both species. However, it is well known that the C-bearing 
molecules freeze out in the central core parts whereas the 
N-bearing molecules survive. Thus, observing the NH3 and 
HC3N molecules we cannot completely reject the possibil- 
ity that the velocity shifts between them are of kinematic 
nature. 

It is clear that the simultaneous observation of the NH3 
(1 ,1) inverse transition and rotational transitions of some 
N-bearing molecules such as, e.g., N 2 H + (1-0) and N 2 D + 
(1-0) would give a more stringent test. The main obstacle 
on this way is that the laboratory frequencies of N 2 H + and 
N 2 D + are known with accuracy not better than 14 m s _1 
(Paper I, Sect. 4. 2). Using the N 2 H + rest frequency from 
the Cologne Database for Molecular Spectroscopy (CDMS; 
Miillcr et al. 2005) and observing with the Nobeyama 45- 
m telescope we obtained a velocity shift between N 2 H + 
and NH 3 of 23.0 ± 3.4 m s" 1 in L1498, 24.5 ± 4.3 m s" 1 
in L1512, and 21.0 ± 5.1 m s" 1 in L1517B (Paper I). We 
note that similar shifts are indicated for the central parts 
of L1498 and L1517B in Fig. 11 from Tafalla et al. (2004) 
who observed these targets with the 30-m IRAM telescope 
and used the N 2 H+ rest frequency very close to the CDMS 
value. 

In the prestellar core L183, Pagani et al. (2009) observed 
N 2 H+ (1-0) and N 2 D+ (1-0) and found that along two 
perpendicular cuts the two lines are systematically shifted 
with respect to NH3 (1,1). The velocity offset is about 28 
m s _1 , when adopting the CDMS catalogue rest frequencies 
(Paper I; Molaro et al. 2009), but AV « 40 m s" 1 if the 
frequency of N 2 H + is taken from Caselli et al. (1995). 

Thus, almost the same velocity offset between the 
NH3 inversion transition and rotational lines from differ- 
ent molecular species (HC3N, N 2 H + , N 2 D + ), different tar- 
gets, and different telescopes is observed. Tentatively, this 
supports our initial hypothesis that the equivalence prin- 
ciple of local position invariance may break in low-density 
environment. Obviously, for more definite conclusions new 
laboratory measurements of the rest frequencies and new 
observations involving other targets and other rotational 
transitions are required. 

In this context we note that the matter density in dif- 
fuse interstellar clouds, where the UV lines of the molecu- 
lar hydrogen H 2 are observed, is about 5 orders of mag- 
nitude lower than in molecular cores. This implies that 
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the expected value of A/i//i based on the H2 line posi- 
tion measurements should not be less than 3 x 1CP 8 . This 
makes future experiments with CODEX and ESPRESSO 
spectrograph for extremely stable Doppler measurements 
(Pasquini et al. 2009) of a great importance since cur- 
rent limits on A/i//i obtained with the UVES/VLT and 
HERES/Keck are at the level of a few ppm (e.g., Malec et 
al. 2010). 

5.2. Frequency shifts caused by external fields 

In this section we consider possible systematic effects 
caused by external fields: we start with a static electric 
field, then consider a thermal background radiation, and, 
finally, discuss a static magnetic field. We deduce all esti- 
mates for the inversion transition in ammonia. Rotational 
transitions can be treated in a similar way. 

5.2.1. Stark effect induced by static electric fields 

Static electromagnetic fields are always present and can 
affect the observed frequencies of microwave transitions in 
the laboratory and in astrophysical environments (see, e.g., 
Bethlem et al. 2008). 

The static Stark displacement of the inversion frequency 
in NH3 is given by (Townes & Schawlow 1955) 



Suiiriv — 2 



(V) 2 T 2 



= 2 



MK 



1 2 



J(J+ 1) 



V 2 T 2 



(10) 



Here fi = h/2ir is the reduced Planck constant, (V) is the 
electric dipole amplitude for the NH3 (1,1) inversion transi- 
tion, P = 1.42 D is the molecular dipole moment, T is the 
electric field strength, J, M, and K are, respectively, the 
molecular angular momentum and its projections on the 
external electric field and on the molecular axis. Averaging 
over the projection of M, one findsB 



SV S t 8wh 
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where <5Vst is the line-of-sight velocity offset caused by the 
Stark effect. The largest offset occurs for the K — J tran- 
sitions. In this case the shift of 5Vst = 1 m s _1 requires a 
field of T « 2 V cm -1 . For laboratory conditions, this is a 
sufficiently large field which is easily controlled. For dense 
interstellar molecular clouds, such a field would accelerate 
electrons to ultra-relativistic energies (~ 10 10 eV) within 
their mean free path (4 ~ 5 x 10 14 /n H ., cm). Thus, we con- 
clude that static electric fields cannot cause any systematic 
effects on this scale. 

We note that the collision shift when molecules get close 
to each other and the NH3 (1,1) inversion line is shifted due 
to the electric dipole-induced interactions between collision 
partners is also less than 1 m s _1 for the physical condi- 
tions realized in the most accurate laboratory studies with 
molecular fountains and in cold interstellar molecular cores 
(Bethlem et al. 2008). 



5 The projection of M ranges between — J and J. Assuming 
equal population of these levels (all of them have the same en- 
ergy), we have (M 2 ) = ^ Y?m=-j M * = J ( J + 1)A 



5.2.2. Stark effect induced by black body radiation 

Now we estimate the frequency shift due to the black 
body radiation (BBR) at a given radiation temperature 
T. According to Farley & Wing (1981), the BBR-induced 
Stark shift is given by: 



c 3 h 4 
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where k is the Boltzmann constant, and c is the speed of 
light. In Eq. (ri2"|) . we take into account that the levels of 
the inversion doublet are shifted from each other and sum 
over Cartesian components of the BBR field. The universal 
function d> has the form 
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(13) 



and is restricted to the interval — 2 < (f> < 3. At the radi- 
ation temperature T > 3 K, we are interested in the limit 
y = (huJi nv /kT) <C 1, where 



(14) 



Using this expression and averaging again over the projec- 
tion of M, we obtain: 
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We see, that the BBR shift grows with temperature and 
it is larger for a warm laboratory environment as compared 
with the interstellar 3 K background field. However, even 
for T ~ 300 K the BBR shift for the inversion transition is 
negligible, Suj inv /uj inv ~ 1.5 x 10 -13 . The BBR shift for the 
J = 2 — 1 rotational transition of the HC3N molecule is of 
the same order of magnitude. 

5.2.3. Zeeman effect induced by static magnetic fields 

Finally, we estimate the effect of a static magnetic field. 
Dense molecular cores have weak magnetic fields with typ- 
ical strengths of B ~ 10 /iG or less (Crutcher et al. 2010). 
For example, CCS Zeeman observations of L1498 set an up- 
per limit on B = 48 ±31 /xG for the line of sight component 
of the magnetic field (Levin et al. 2001). The line shapes 
can be sensitive to random motions at the Alfven velocity 
associated with the magnetic field B, v A = B/y/Airp (here 
B is the strength in G of the line-of-sight component of 
the magnetic vector, p is the ion density in g cm -3 , and 
v A is in cm s _1 ). The presence of the magnetic field can be 
felt by the neutral material through collisions with charged 
particles. If we assume that in our observations the non- 
thermal line widths are due to Alfvenic turbulent motions 
whose kinetic energy is equal to the static magnetic field 
energy, so that v A = v^wturb, then an order of magnitude 
estimate of the magnetic field strength at n H2 ~ 10 5 cm~ 3 , 
and Vturb < 90 m s" 1 (see Sect. [3]) provides an upper limit 
on B < 20 /iG. Since in this estimate we substitute the ion 
density by the total gas density, the strength of the static 
magnetic field in the interiors of LI 498, L1512, L1517B, 
and L1400K should be considerably lower than this upper 
limit. This is more than five orders of magnitude smaller 
than the unscreened Earth's magnetic field which is about 
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Table 6. Mean velocity offsets AV (in m s *) between the HC3N (2-1) and NH 3 (1,1) lines. The values in parentheses 
arc the la uncertainties. 





L1498 


L1512 


L1517BC 


L1400K 


Average 




<AV.) 


(AV a ) 


(AV.) 


(AV a ) 


(AV) 


(AV a ) 


(AV) 


(AVa) 


(AV) 


(AV a ) 


Unweighted 


23.6(2.9) 


22.7(3.6) 


27.0(2.8) 


26.6(2.0) 


47.9(3.1) 


46.9(2.7) 


9.3(4.3) 


8.0(3.4) 


27.9(2.6) 


27.1(2.5) 


Weighted 


24.8(2.2) 


24.8(2.4) 


28.5(1.9) 


27.2(1.8) 


49.2(3.0) 


48.4(2.7) 


9.8(4.3) 


8.5(3.4) 


29.6(2.2) 


28.6(2.1) 


M-estimate 


26.0(1.6) 


25.8(1.7) 


27.8(2.3) 


28.0(1.8) 


47.9(3.9) 


46.9(3.3) 


9.8(4.5) 


8.5(3.4) 


27.5(2.4) 


27.3(1.5) 


Median 


24.8 


24.3 


26.5 


26.5 


48.3 


46.8 


10.5 


9.0 


26.5 


26.5 


Sample size 


8 


8 


13 


13 


6 


6 


4 


4 


31 


31 



half a Gauss at the surface of the Earth. Typical laboratory 
magnetic shields reduce the Earth's field by 2-3 orders of 
magnitude to a mG level. 

The order of magnitude estimate of the Zeeman split- 
ting for NH3, where magnetic interaction is dominated by 
3 protons since the nitrogen (/-factor is small (w 0.4), can 
be obtained as follows. For a single proton, we have 

5uj inv = /j, n g p I z B , (16) 

where fi n is the nuclear magneton, g p = 5.6 is the proton 
g-factor, I z is the proton spin and B is the strength of 
the magnetic field. Roughly, the maximum shift is realized 
when all three proton spins are aligned (ortho-NH^), so we 
can use the above equation and put I z — I z ,i + I z .2 + I z .3 = 
3/2. With fj, n = 5.05 x 10~ 24 erg G" 1 and B = 10~ 3 G, we 
obtain 5u- mv w 6 Hz. Since the (J,K) = (1, 1) line belongs 
to the para-species of the ammonia molecule, where the 
proton spins are not aligned, the expected frequency shift 
is even lower for the transition observed by us. 

One should also keep in mind that the Zeeman shifts 
turn to zero after averaging over magnetic quantum number 
M . Therefore, in the first approximation, magnetic fields 
induce a line broadening, but not a frequency shift. Thus, 
the systematic shift of the ammonia inversion transition 
of about <~ 2 kHz detected in the starless molecular cores 
cannot be explained by a static magnetic field effect. 

We conclude that static fields and thermal background 
radiation produce systematic velocity shifts much smaller 
than 1 m s _1 . Thus, field effects do not contribute signifi- 
cantly to the systematic error of our observations. 

6. Conclusions 

The present paper continues our studies where we use the 
ammonia inversion transition in conjunction with low-lying 
rotational transitions of other molecules to probe the de- 
pendence of the electron-to-proton mass ratio /j, on the 
ambient matter density. As experimental tests, we sug- 
gest to observe dense prestellar molecular cores located in 
the disk of the Milky Way. Here we report on mapping of 
four molecular cores (L1498, L1512, L1517B, and L1400K) 
in two molecular transitions of NH3 (J,K) = (1,1) and 
HC3N J = 2 — 1 with high-spectral resolution at the 100-m 
Effclsberg telescope. The main results are as follows. 

1. The completely resolved hfs components of the NH3 
(1,1) and HC3N (2-1) transitions allow us to determine 
the line centers with the precision of ~ 1 m s _1 , which 
is comparable with laboratory uncertainties of the rest 
frequencies of these transitions. 



2. In two cores, L1498 and L1512, with lowest Dopplcr 
noise, we obtain a statistically significant positive ve- 
locity offset between the rotational HC3N and inversion 
NH3 transitions of 26.9 ± 1.2 sta t ± 3.0 sys m s _1 . In two 
other cores, L1517B and L1400K, the velocity offsets 
are, respectively, higher and lower, but these deflections 
can be explained by the observed kinematic structure 
in L1517B and L1400K. If we assume that the mea- 
sured velocity offset is caused by the electron-to-proton 
mass ratio variation, then A/i/fi — 26 ± l sta t ± 3 sys 
ppb. The non-zero A/u implies that at deep interstel- 
lar vacuum the electron-to-proton mass ratio increases 
by ~ 3 x 10~ 8 as compared with its terrestrial value. 

3. The reproducibility of the velocity offset at the same 
facility (Effclsberg telescope) on the year-to-year base 
is very good, except for the L1400K target where the 
central point is probably an outlier. 

4. We show that the effects of static electric and magnetic 
fields and the BBR-induced Stark shifts are not larger 
than 1 m s _1 , and can be neglected in the total error 
budget. 

5. The results obtained tentatively support the hypothesis 
that in low-density environment the equivalence prin- 
ciple may break. This may be a consequence of the 
chameleon-like scalar field. New laboratory measure- 
ments of the rest frequencies and new observations in- 
volving other rotational transitions and other targets 
(e.g., Levshakov et al. 2010b) are required to reach more 
definite conclusions. 
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